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Abstract 

The basic theory of Differential Galois and in particular Morales- 
Ramis theory is reviewed with focus in analyzing the non-integrability 
of various problems of few bodies in Celestial Mechanics. The main 
theoretical tools are: Morales-Ramis theorem, the algebrization me- 
thod of Acosta-Blazquez and Kovacic's algorithm. Morales-Ramis 
states that if Hamiltonian system has an additional meromorphic in- 
tegral in involution in a neighborhood of a specific solution, then the 
differential Galois group of the normal variational equations is abelian. 
The algebrization method permits under general conditions to recast 
the variational equation in a form suitable for its analysis by means of 
Kovacic's algorithm. We apply these tools to various examples of few 
body problems in Celestial Mechanics: (a) the elliptic restricted three 
body in the plane with collision of the primaries; (b) a general Hamil- 
tonian system of two degrees of freedom with homogeneous potential 
of degree —1; here we perform McGehee's blow up and obtain the 
normal variational equation in the form of an hyper geometric equa- 
tion. We recover Yoshida's criterion for non-integrability. Then we 
contrast two methods to compute the Galois group: the well known, 
based in the Schwartz-Kimura table, and the lesser based in Kovacic's 
algorithm. We apply these methodology to three problems: the rect- 
angular four body problem, the anisotropic Kepler problem and two 
uncoupled Kepler problems in the line; the last two depend on a mass 
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parameter, but while in the anisotropic problem it is integrable for 
only two values of the parameter, the two uncoupled Kepler problems 
is completely integrable for all values of the masses. 

Keywords and Phrases. Algebrization method, Differential Galois 
Theory, Celestial mechanics, Kovacic's algorithm, Kimura's theorem, 
Morales-Ramis theory, n— body problems, non-integrability. 
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1 Introduction 

In this paper we analyze the integrability of some Hamiltonian systems of 
two degrees of freedom related with few body problems. This can be made 
through the analysis of the linearization of the Hamiltonian system, that 
is, variational equations and normal variational equations. In 1982 Ziglin 
(|25j) proved a non-integrability theorem using the constraints imposed on 
the monodromy group of the normal variational equations along some integral 
curve by the existence of some first integrals. This is a result about branching 
of solutions: the monodromy group express the ramification of the solutions 
of the normal variational equation in the complex domain. 

We consider a complex analytic symplectic manifold M of dimension 2n 
and a holomorphic hamiltonian system Xh denned over it. Let V be the 
Riemann surface corresponding to an integral curve z = z(t) (which is not an 
equilibrium point) of the vector field Xh- Then we can write the variational 
equations (VE) along V, 



Using the linear first integral dH(z(t)) of the VE it is possible to reduce this 
variational equation (i.e. to rule out one degree of freedom) and to obtain 
the so called normal variational equation (NVE) that, in some adequate 
coordinates, we can write, 
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where, as usual, 

J=( ° 1 
J {-10 

is the square matrix of the symplectic form. (Its dimension is 2(n — 1)). 

In general if, including the hamiltonian, there are k analytical first inte- 
grals independent over T and in involution, then, in a similar way, we can 
reduce the number of degrees of freedom of the VE by k. The resulting 
equation, which admits n — k degrees of freedom, is also called the normal 
variational equation (NVE). Then we have the following result (|25j). 



Theorem (Ziglin). Suppose that the hamiltonian system admits n — k ad- 
ditional analytical first integrals, independent over a neighborhood of V (but 
not necessarily on V itself) We assume moreover that the monodromy group 
of the NVE contain a non-resonant transformation g. Then, any other el- 
ement of the monodromy group of the NVE send eigendirections of g into 
eigendirections of g. 

We recall that a linear transformation g G Sp(m, C) (the monodromy 
group is contained in the symplectic group) is resonant if there exists integers 
ri, ...,r m such that A^ 1 • • • A r ™ = 1 (where we denoted by \ the eigenvalues 
of <?)■ 

Later, Morales and Ramis in 2001 improved the Ziglin's result by means 
of differential Galois theory (see [H] and see also [H]), arising in this way 
the so-called Morales-Ramis Theory. This theory will be explained in section 
[3] of this paper. 

There are a lot of papers and books devoted to analyze three body prob- 
lems (see [T7] and references therein). Therefore, an special kind of three 
body problem is the so-called Sitnikov problem, which has been deeply an- 
alyzed using Morales-Ramis theory in [3j [HI EH]- Another cases of three 
body problems has been studied, also by means of Morales-Ramis theory, in 
[6l [7J. There are a lot cases in which the variational equation falls in Rie- 
mann differential equation or hypergeometric differential equation. In this 
cases has been used satisfactory the Kimura-Schwartz table (see [10J ) , which 
was improved by Morales in [XT] . 

In this paper, we analyze the non-integrability of some celestial mechanics 
problems such as the collinear restricted elliptic three-body problem, rect- 
angular 4 body problem and the anisotropic Kepler problem. The approach 
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used here is by means of Morales Ramis theory contrasting the Kovacic's 
algorithm with Kimura's theorem, but obtaining the same results. 



2 Differential Galois Theory 

Our theoretical framework consists of a well-established crossroads of Dy- 
namical Systems theory, Algebraic Geometry and Differential Algebra. See 
[TT] or [18] for further information and details. Given a linear differential 
system with coefficients in C(t), 

z = A{t)z, (1) 

a differential field L D C(t) exists, unique up to C(t)-isomorphism, which 
contains all entries of a fundamental matrix \1/ = [tpx, . . . , ip n ] of ([1]). More- 
over, the group of differential automorphisms of this field extension, called 
the differential Galois group of ([I]), is an algebraic group G acting over the C- 
vector space (^i, . . . , ip n ) of solutions of ([1]) and containing the monodromy 
group of ([1]). 

It is worth recalling that the integrability of a linear system (JT|) is equiva- 
lent to the solvability of the identity component G° of the differential Galois 
group G of ([T])- in other words, equivalent to the virtual solvability of G. 

It is well established (e.g. [3j [12]) that any linear differential equation 
system with coefficients in a differential field K 



d_ ( & \ / a{t) b(t) \ ( 6 
dt\&) V <M <*) J U 



(2) 



by means of an elimination process, is equivalent to the second-order equation 

I - (<t) + d(t) + M b _ ^(t) + 6(t)c(*) - a(t)d(t) - f = 0, (3) 

where ^ := ^x- Furthermore, any equation of the form z — 2pz — qz = 0, can be 
transformed, through the change of variables z = ye^ p , into ij = —ry, r satisfying 
the Riccati equation p = r + q+p 2 . This change is useful since it restricts the study 
of the Galois group of y = —ry to that of the algebraic subgroups of SL(2, C). This 
last procedure will be used later in section algorithmic approach. 

A natural question which now arises is to determine what happens if the coeffi- 
cients of the differential equation are not all rational. A new method was developed 
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in [JJ, in order to transform a linear differential equation of the form x = r(t)x, 
with transcendental or algebraic non-rational coefficients, into its algebraic form 
- that is, into a differential equation with rational coefficients. This is called the 
algebrization method and is based on the concept of Hamiltonian change of vari- 
ables [4J. Such a change is derived from the solution of a one-degree-of- freedom 
classical Hamiltonian. 

Definition 1 (Hamiltonian change of variables). A change of variables r = r(i) 
is called Hamiltonian if (r(f), f(t)) is a solution curve of the autonomous Hamil- 
tonian system Xh with Hamiltonian function 



Theorem 1 (Acosta-Blazquez algebrization method [4J). Equation x = r{t)x is 
algebrizable by means of a Hamiltonian change of variables r = r(t) if, and only 



The next intended step, once a differential equation has been algebrized, is 
studying its Galois group and, as a causal consequence, its integrability. Concern- 
ing the latter, and in virtue of the invariance of the identity component of the 
Galois group by finite branched coverings of the independent variable (Morales- 
Ruiz and Ramis, [141 Theorem 5]), it was proven in [4J Proposition 1] that the 
identity component of the Galois group is preserved in the algebrization mecha- 
nism. 

The final step is analyzing the behavior of t = oo (or r = oo) by studying 
the behavior of r) = through the change of variables n = 1/t (or r) = 1/t) 
in the transformed differential equation, i.e. t = oo (or r = oo) is an ordinary 
point (resp. a regular singular point, an irregular singular point) of the original 
differential equation if, and only if, rj = is one such point for the transformed 
differential equation. 



H = H(r,p) = y + V(t), for some V G C(r). 




Furthermore, the algebraic form of x = r(t)x is 




□ 
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3 Morales- Ramis Theory 



Everything is considered in the complex analytical setting from now on. The heuris- 
tics of the titular theory rest on the following general principle: if we assume system 



"integrable" in some reasonable sense, then the corresponding variational equa- 
tions along any integral curve r = {z(t) : t £ /} of ([5|), defined in the usual man- 
ner 



must be also integrable - in the Galoisian sense of the last paragraph in [2j We 
assume T, a Riemann surface, may be locally parameterized in a disc / of the 
complex plane; we may now complete T to a new Riemann surface T, as detailed 
in [HI §2.1] (see also [HI §2.3]), by adding equilibrium points, singularities of the 
vector field and possible points at infinity. Linearization defines a linear connection 
over T called the variational connection VEp- and Gal (VEp) is its Galois differential 
group which contains the Zariski closure of the monodromy group Mon (VEp) . In 
practice the normal variational equations NVE^ = TM\r/TT are analyzed, the 
variational equation along the solution being reducible. 

The aforementioned "reasonable" sense in which to define integrability if sys- 
tem ([5]) is Hamiltonian is obviously the one given by the Liouville- Arnold Theorem 
(see [HO [22]), and thus the above general principle does have an implementation: 

Theorem 2 (J. Morales- Ruiz & J. -P. Ramis, 2001). Let H be an n- degree- of- 
freedom Hamiltonian having n independent rational or meromorphic first integrals 
in pairwise involution, defined on a neighborhood of an integral curve T. Then, the 
identity component Gal (VEp) is an abelian group (i.e. Gal (VEp) is virtually 



The disjunctive between meromorphic and rational Hamiltonian integrability 
in Theorem [2] is related to the status of t = oo as a singularity for the normal 
variational equations. More specifically, and besides the non-abelian character of 
the identity component of the Galois group, in order to obtain Galoisian obstruc- 
tions to the meromorphic integrability of H the point at infinity must be a regular 
singular point of flVErD (for example Hypergeometric and Riemann differential 
equations). On the other hand, for there to be an obstruction to complete sets of 
rational first integrals, t = oo must be a irregular singular point. See [141 Corollary 
8] or [11, Theorem 4.1] for a precise statement and a proof. 




(5) 



i = X'(z(t))t, 



(VE r ) 



abelianj. 
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Different notions of integrability correspond to classes of admissible first in- 
tegrals, for instance rational, meromorphic, algebraic, smooth, etc. For non- 
integrability within the real-analytic realm, a popular testing is Melnikov integral 
whose isolated zeros give transversal homoclinic intersections and in under proper 
hypothesis, chaos. Galoisian obstruction to integrability based in Morales-Ramis 
theory has shown to be equivalent to the presence of isolated zeros of Melnikov in- 
tegral, for a class of Hamiltonian systems with two degrees of freedom with saddle 
centers (see [131 [23]). Also high order variational equations have been studied in 
this context (see [16]). 

Remark 1. In order to analyze normal variational equations, a standard procedure 
is using Maple, and especially commands dsolve and kovacicsols. Whenever 
the command kovacicsols yields an output "[ ]", it means that the second-order 
linear differential equation being considered has no Liouvillian solutions, and thus 
its Galois group is virtually non-solvable. For equations of the form y = ry with 
r £ C(x) the only virtually non-solvable group is SL(2,C). In some cases, more- 
over, dsolve makes it possible to obtain the solutions in terms of special functions 
such as Airy functions, Bessel functions and hypergeometric functions, among 
others (J^)- There is a number of second-order linear equations whose coefficients 
are not rational, and whose solutions Maple cannot find by means of the com- 
mands dsolve and kovacicsols alone; this problem, in some cases, can be solved 
by the stated algebrization procedure. Another difficulty is when appears parame- 
ters in the differential equation, then almost always kovacicsols wrong, for this 
reason we present in following section the Kovacic 's algorithm to be used later. 

4 Kovacic's Algorithm 

This algorithm is devoted to solve the RLDE (reduced linear differential equation) 
£" = r£ and is based on the algebraic subgroups of SL(2,C). For more details see 
[9]. Improvements for this algorithm are given in |19j . where it is not necessary 
to reduce the equation. Another improvement is given in [8], which is a compact 
version to implement in computer systems. Here, we follow the original version 
given by Kovacic in [9], which is the same version given in [3]. 

Theorem 3. LetG be an algebraic subgroup o/SL(2,C). Then one of the following 
four cases can occur. 

1. G is triangularizable. 

2. G is conjugate to a subgroup of infinite dihedral group (also called meta- 
abelian group) and case 1 does not hold. 
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3. Up to conjugation G is one of the following finite groups: Tetrahedral group, 
Octahedral group or Icosahedral group, and cases 1 and 2 do not hold. 

I G = SL(2,C). 

Each case in Kovacic's algorithm is related with each one of the algebraic 
subgroups of SL(2, C) and the associated Riccatti equation 

e' = r -e 2 = (yfi-o) (VF + 0) , e = j. 

According to Theorem [3j there are four cases in Kovacic's algorithm. Only for 
cases 1, 2 and 3 we can solve the differential equation the RLDE, but for the case 
4 we have not Liouvillian solutions for the RLDE. It is possible that Kovacic's 
algorithm can provide us only one solution so that we can obtain the second 
solution (£ 2 ) through 

«•-*/!• (6) 

Notations. For the RLDE given by 

-^2= r£, r=-, 8,t€C[x], 

we use the following notations. 

1. Denote by T' be the set of (finite) poles of r, T' = {c G C : t(c) = 0}. 

2. Denote by T = T' U {oo}. 

3. By the order of r at c S T', °{r c ), we mean the multiplicity of c as a pole of 
r. 

4. By the order of r at oo, o (r^) , we mean the order of oo as a zero of r. That 
is o (r-oo) = deg(t) - deg(s). 

4.1 The four cases 

Case 1. In this case [y/r\ c and [y^oo means the Laurent series of y/r at c and 
the Laurent series of \fr at oo respectively. Furthermore, we define e(p) as follows: 
if p € T, then e (p) G {+,—}• Finally, the complex numbers a+,a~,a+ , a~ will 
be defined in the first step. If the differential equation has not poles it only can 
fall in this case. 

Step 1. Search for each c G T' and for oo the corresponding situation as 
follows: 
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(co) If o (r c ) = 0, then 
(ci) If o (r c ) = 1, then 
(C2) If o (r c ) = 2, and 



[Vr] c = 0, o£ = 0. 

[Vr] c = o, «± = i. 

r = • • • + 6(x — c)~ 2 + • • • , then 



W r \c = °> a c = " 



2 

(c 3 ) If o (r c ) = 2v > 4, and 

r = (a (x - c)~ v + ... + d (x - c)" 2 ) 2 + b(x - c)- (v+1) + ■■■ , then 

[VF] C = a(x - c)~ v + ... + d(x - c)- 2 , o£ = ^ (±^ + u ) • 

(001) If o (roo) > 2, then 

(002) If o (roo) = 2, and r = ■ ■ ■ + 6x 2 + • • • , then 

Lv r Joo = ' °°o = g • 

(003) If o (r^) = — 2f < 0, and 

r = (ax 11 + ... + d) 2 + bx v ~ 1 H , then 

[Vr] OQ = ax v + ... + d, and a±=^(±^-vj. 

Step 2. Find Z? / defined by 

D = \m € Z + : m = a*°> - V (e (p)) peT \ . 

I cgT' J 

If D = 0, then we should start with the case 2. Now, if #D > 0, then for each 
m e D we search oj £ C(x) such that 

" = e (00) [v^ + £ (e (c) [Vr] c + <&\x - c)' 1 ) . 
cex' 
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Step 3. For each m G D, search for a monic polynomial P m of degree m with 

Pi + 2uP' m + (u/ + cu 2 - r)P m = 0. 

If success is achieved then £1 = P m e§ ^ is a solution of the differential equation 
the RLDE. Else, Case 1 cannot hold. 

Case 2. Search for each c£T' and for oo the corresponding situation as 
follows: 

Step 1. Search for each c G T' and oo the sets E c ^% and -Eqo 7^ 0. For each 
c e T' and for 00 we define E c d7L and C Z as follows: 

( Cl ) If o (r c ) = 1, then E c = {4} 

(02) If o (r c ) = 2, and r = • • • + 6(x — c)~ 2 + • • • , then 

£ c = {2 + kVTT4b : fc = 0, ±2} . 

(c 3 ) If o (r c ) = i; > 2, then £ c = {v} 

(001) If o ( roo ) > 2, then = {0, 2, 4} 

(002) If (^oo) = 2, and r = ■ ■ ■ + 6x 2 + • • • , then 

^ = {2 + kVl + 46 : fc = 0, ±2} . 

(003) If o (r^) = v < 2, then = {v} 
Step 2. Find Z? / defined by 



D = < m £ Z + : m = 



If D = 0, then we should start the case 3. Now, if j^D > 0, then for each m € D 
we search a rational function 9 defined by 
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Step 3. For each m £ D, search a monic polynomial P m of degree m, such 
that 

P'Z + 39P^ + (39' + 36 2 - Ar)P' m + (6" + 3##' + 6 3 - Ar9 - 2r') P m = 0. 

If P m does not exist, then Case 2 cannot hold. If such a polynomial is found, set 
4> = + P' IP and let oj be a solution of 



+ </>w + -(>' + 2 - 2r) = 0. 
Then = is a solution of the differential equation the RLDE. 



Case 3. Search for each c£ T' and for oo the corresponding situation as 
follows: 

Step 1. Search for each c G T' and oo the sets E c ^§ and / 0. For each 
c G T' and for oo we define £7 C C Z and £?oo C Z as follows: 

(ci) If o(r c ) = l, then ^ = {12} 

(C2) If o (r c ) = 2, and r = h 6(x - c)~ 2 H , then 

£ c = 1 6 + kVl + 46 : fc = 0, ±1, ±2, ±3, ±4, ±5, ±6} . 

(oo) If o (roo) = v > 2, and r = • • • + te 2 + • • • , then 

f lot 1 

Eoo = ie + —VT+4b : jfe = 0,±l,±2,±3,±4,±5,±6L ne{4,6,12}. 



Step 2. Find Z? 7^ defined by 



00 



:ST' / J 



In this case we start with n = 4 to obtain the solution, afterwards n = 6 and 
finally ra = 12. If D = 0, then the differential equation has not Liouvillian solution 
because it falls in the case 4. Now, if #D > 0, then for each m G D with its 
respective n, search a rational function 
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ceT' 
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and a polynomial S denned as 

s=n(x-c). 

cex' 

Step 3. Search for each m £ D, with its respective n, a monic polynomial 
P m = P of degree m, such that its coefficients can be determined recursively by 

P_l = 0, P n = -P, 

Pi-X = -SP( - ((n - i) 5' - 50) Pi - (n - i) (t + 1) 5 2 rP 4+1 , 

where i E {0, 1 . . . , n — 1, n}. If P does not exist, then the differential equation has 
not Liouvillian solution because it falls in Case 4. Now, if P exists search u such 
that 

^(n-*)! 

then a solution of the differential equation the RLDE is given by 
where w is solution of the previous polynomial of degree n. 



4.2 Some remarks on Kovacic's algorithm 

Along this section we assume that the RLDE falls only in one of the four cases. 

Remark 2 (Case 1). If the RLDE falls in case 1, then its Galois group is given 
by one of the following groups: 

11 e when the algorithm provides two rational solutions or only one rational solu- 

tion and the second solution obtained by ^ has not logarithmic term. 

1 
1 

this group is connected and abelian. 

12 G/c when the algorithm provides only one algebraic solution £ such that £ fc 6 

C(x) andi k ~ l C(z). 

1 J : X is a k-root of the unity, d G C 



A 

this group is disconnected and its identity component is abelian. 
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13 C* when the algorithm provides two non-algebraic solutions. 

C * = {(o c-) :CeC 
this group is connected and abelian. 

14 C + when the algorithm provides one rational solution and the second solution 

is not algebraic. 



c+ = {(o l) :d€C }' £ eC (*)' 



this group is connected and abelian. 

15 C* x C + when the algorithm only provides one solution £ such that £ and its 

square are not rational functions. 

C*xC+ = |^ J 1 ^ :c€C*,d€c|, C£C(x), f<£C(x). 

This group is connected and non- abelian. 

16 SL(2, C) if the algorithm does not provide any solution. This group is connected 

and non-abelian. 

Remark 3 (Case 2). If the RLDE falls in case 2, then Kovacic's Algorithm can 
provide us one or two solutions. This depends on r as follows: 

III if r is given by 

2(j)' + 2(f) - (j) 2 



r 

4 

then there exist only one solution, 

112 if r is given by 

20' + 20 - (j? 

T + A ' 

then there exists two solutions. 

113 The identity component of the Galois group for this case is abelian. 

Remark 4 (Case 3). If the RLDE falls in case 3, then its Galois group is given 
by one of the following groups: 
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III1 Tetrahedral group when to is obtained with n = 4. This group of order 24 
is generated by 



e~ \ 1 ( n hjH \ (\ 1 \ 

' 3l 2e3 - 1 )[2 k£Z - 



kiri 

e" — 



1112 Octahedral group when to is obtained with n = 6. This group of order 48 
is generated by 

e-^J' 2 £4 ( C2 * eZ - 

1113 Icosahedral group w/ien w is obtained with n = 12. This group of order 
120 is generated by 



e~ \ /0V 
e - ^ i ' \V -0 



, fc G Z, 



&em<? and V defined as 

1 / 3fc7ri 2kiri kiri \ 1 / Zkni 2kTvi kiri \ 

4> = - \^e~ - e~ + 4e~ - 2J , V = g (^e~ + 3e~ - 2e~ + lj 
III4 TTte identity component of the Galois group for this case is abelian. 



5 Applications 

5.1 The collinear restricted elliptic three-body prob- 
lem 

Let two primaries of mass 1 — fi and fi move along the x-axis, its positions being 
X\ = —fir, X2 = (l — fi)r where fi = m<il{vn,\-\-m<2). Suppose the primaries perform 
an elliptic collision motion 

r = 1 — cos E 
t = E-sinE 

where E is the elliptic anomaly and we choose units of time and length such that 
the maximum distance between the primaries is unit and the mean motion is one. 
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The equations of motion of a massless particle in a fixed plane containing the line 
of the primaries is 



(l- / u)(x-xi) n{x — xi) 



\x — 

- fi)(x+ fir) 
\x + //r| 3 



\x — 
fi(x 



X2 



(l-fi)r) 



(l-/i)r| 



(7) 



where its position is x G C. System (|7|) is two degrees of freedom time-dependent 
Hamiltonian system. Some general results are known for time-dependent, one de- 
gree of freedom (see for example [3]). For the present we take an ad- hoc procedure: 
We will perform several changes of variables in order to obtain the desired form of 
equations of motion. Firstly, perform a change to pulsating coordinates 



then (El) transforms into 



Using the elliptic anomaly E as independent variable, 

d 1 d 



1 + A*T 



(8) 



0) 



dt r dE 



yields 



1 d 

r rdE 
d 

dE 



1 de 

r d£ 
1 d£ 

7Ze 



1 dr 1 d£ 
+ 2 - 

r di? r d.E 

2 sin 75 d£ 



d£ 



ie+Mi 



+ 



where the potential function is 

m = he +^7, + iH i T7T- ( 10 ) 

Developing the left hand side of the previous ode we obtain 
sin E d£ 1 d 2 £ 2sin£ d£ 1 

ld 2 £ sinEd^ 1 
rd£ 2 + r 2 dE " r 2 Vi '^ 
d 2 £ sinEde 1 
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In summary, 

d 2 £ sinE d£ 1 



The equation (jlip can be analytically extended to the whole complex E'-plane 
except for singularities at the point on the real axis E = ±7r, ±2ir, . . ., and also 
has singularities due to collisions with the binaries £ = — fj,, 1 — fj,. 

The critical points of (|1U|) are given by the classical Eulerian and Lagrangian 
points satisfying Vfi^^J = 0, i = 1, 2 . . . , 5. Let B denote the Hessian 

o _ ( ^11 ^12 
V ^21 ^22 

evaluated at any of the points Lj. The linearization of ljlip at Lj is 

d 2 j sinE d£ _ Btj (u) 
dE 2 1-cosEdE 1-cosfT 

The above procedure can be seen as the linearization of the lifted system 



1, 



dE 
ds 

ds 

dv sinE di B£ 



ds 1 — cos E dE 1 — cos E 

where E is considered mod 2tt, along any of the periodic orbits £ = ^n, i = 1,2, 3, 
E = s. 

The following properties of matrix B are well known (for details see [20]): For 
collinear configurations Li,L2,L^, B = diag(/«i, K2), with k\ > and K2 < for 
all values of the mass parameter [i, the exact values depend on the root of Euler's 
quintic equations. 

In this case, the variational equations (|12p split 



d 2 £i sinE d£j , . , 

77^+1 ?^l = r^^' where.? = 1,2. 13 

a£/ z 1 — cos E dE 1 — cos E 

Theorem 4. For the collinear elliptic restricted three body problem in the plane 
fH let 



2 1 ie+Mi ie-i+Mi 

and Zet J i/ie se£ 0/ exceptional mass parameters [i £ (0, 1) suc/i i/iaf: 
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(i) £*(//), * = 1,2,3 belongs to any of the three solution curves of the equation 
defining the collinear configurations 

^(^(M),0) = 0; 

-oo<£l(n)<-n, -n<£f(ji) <l-n, 1- n < <oo; 

(ii) any of the coefficients 

«J(m) = (£(/*),<>), «5(Ai) = 0(^M,o) 

«5 = ^ (i4) 

where n is an integer. 
Then if \i J , the 'problem is not integrable. 

Proof. The procedure is to algebrize the variational equations (|13p and then apply 
Kovacic's algorithm. We start considering the variational equation, 

a%- sin£ d£j kj . 

K H — = where 7 = 1,2, 

dE 2 1 - cos EdE l-cos£ 

which is transformed in the differential equation 

dE 2 

where 



HE) y(E) (15) 



_ cos^-l+4 Ki 

4(1 — cos-Ej Vl-cosi? 

Now, by theorem [H the equation (|15p is ready to be algebrized. The Hamiltonian 
change of variable is r = r(S) = cosE, where f = — sinS, (f) 2 = sin 2 E = 
1 — cos 2 E so that 

2 i , r - 1 + 4re,- 

q = 1 — r and t = ; r-^. 

4(1 - r) 

The algebraized equation is 

A(T) T dViT) + Af T ~ 1 t! Ki 2 M r)=0, (16) 



dr 2 1-r 2 dr 4(-l + r)(l-r 2 
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and the points 1, — 1 and oo are regular singularities. To apply Kovacic's algorithm, 
see section 14.14 we use the RLDE 

d^ = r{T)rt > r(T) = 4(l-T)2(l+T)2 ' V{T) = W^ 

with Kj 7^ 0, because k\ > and K2 < 0. We can see that T = {— 1,1, oo} and 
that the equation (|17p could fall in any of four cases of Kovacic's algorithm, now 
expanding r(r) in partial fractions we have that 

. . 8k, — 3 4k,- — 3 3 4k,- — 3 



16(1 -t) 2 16(1 -r) 16(1 + r) 2 16(1 + r) ' 

We start analyzing the case one. The equation (fT7|) satisfy the conditions {c2, ooi}, 
because on = or_i = 2 and or^ = 3, obtaining the expressions 

Wr\-i = Wr\x = [Vr]oo = «£, = 0, a~ = 1, 



3 „ ,- _ l 



± _ 2±- v /8k j +1 



4! ^, Q!^ ^ 

By step two, D = Z + and Kj has the following possibilities: 

Kj = (n + l)(2n + 3), k,- = (n + l)(2n + 1), «j = n(2n+ 1), kj = n(2n - 1) 

which are equivalents to Kj = n(n + l)/2. For each n we can construct u and by 
step three there exists a monic polynomial of degree n in which each solution of 
the differential equation (|17p is given for all n E Z+. 

Following the case two, we expect to find different values of Kj that the pre- 
sented in case one, so that the equation (fT?]) satisfy the conditions {c2,ooi}, be- 
cause on = or_i = 2 and or ro = 3, obtaining the expressions 

Ex = {2, 2 - VT+8kJ, 2 + VT+8k7} , £L X = {1, 2, 3}, £^ = {0,2,4}. 

By step two, D = Z + and Kj is again equivalent to Kj = ra ( n 2 +1 ) , so that we discard 
the case two. 

Finally, following the case 3, we expect to find different values of Kj that the 
presented in case one, but again appear the expression v/T + 8k j , which replaced in 
E c and give us again an equivalent expression to Kj = n ^ n 2 • This means that 
the differential equation (|17p is contained in the Borel group when Kj = w (" 2 +1 ) , 
and it is SL(2, C) when Kj ^ 2 • Therefore, by remark [21 the Galois group is 
virtually abelian for Kj = and unsolvable for Kj ^ n \ n + 1 ) _ □ 
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Sorry, the graphics are not available. 

Figure 1: (a) The curves K,\((i) (up) and (bottom), (b) The curves 

kI((j,) (up) and re|(/i) (bottom). 

An alternative proof based on Kimura's approach is given in the Appendix A.l, 
mainly for contrasting both techniques. 

It is interesting to investigate the exceptional values of the mass parameter 
H £ J such that any of the K%, K2 satisfy the condition (|14|) . Since K2 is negative, 
there are no exceptional values since w ( n 2 +1 ) > for all integers n. Shown in 
Figure Q] are the curves k\ 2 (/■*), for * = 1,2 (k%, K2 are symmetrical with respect 
to fj, = 1/2, so we just consider L\ and L%). For fi = 0, the non-integrability 
test fails since then both k}(0) = 2 • | =3 and k± (0) = 3 • | = 6 (can be verified 
analytically). This is consistent with the fact that for /j, = system is just a Kepler 
problem. The exceptional values \il x satisfying Kiip,^) = 3 • | = 6 and < 0-5 
satisfying «i(//£ 2 ) = 3- 1 = 6 are not satisfied simultaneously for the same value of 
the mass parameter, i.e. / u^ 1 7^ /jll 2 , thus for some of the reference orbits L\ or L 2 , 
the system does not posses an integral in a neighborhood of that orbit, although 
the theorem does not discard the existence of an additional integral locally defined. 



5.2 Homogeneous potential of degree —1 

We consider a general application to a two degrees of freedom simple hamiltonian 
system with homogeneous potential of degree —1 

H = -( P 2 X + P 2 y )-U(x,y), 

We suppose that U{x,y) is defined and is positive for all (x,y) 6 M 2 , except the 
origin. 

The Hamiltonian un polar coordinates becomes 




McGehee's blow up is achieved taking coordinates v = r l ^ 2 p r -, u = r 1 I 2 pq and 
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rescaled time dt = r^/ 2 dr. The equations of motion take the form 

r = rv, 

v' = ^v 2 + u 2 -U(9), 
9' = u, 

u' = --vu + U'(9). 

Where the prime in the left hand side denotes derivatives with respect to r and 
U'(0) denotes derivative with respect to its argument, which causes no confusion. 
System (|18j) leaves invariant the energy surface 

E h = {(r,9,u,v)\r>0, ^(u 2 + v 2 ) = U(9) + rh}, (18) 

which can be extended invariantly up to its boundary, the collision manifold 

A = {(r,9,u,v) | r = 0, ^{u 2 + v 2 ) = U(9)}. (19) 

Since U(9) is periodic, by the mean value theorem, there exists 6 C such that 
U'(6 C ) = 0. Let v c = 2U(6 C ). Then for h < there exists an ejection-collision 
homothetic orbit given explicitly by 9 = 9 C , u = and 

v 2 

r h (r) = -^sech 2 ( Wc r/2), 

Vh(r) = -■u c tanh(v c r/2). 

The variational equations along the homothetic orbit are 

5r' = Vh5r + r^dv, 

5v = Vh5v, 

59' = 5u, 

5v! = --v h 5u + U"(e c )S6 

The last two equations are decouple and constitute the normal variational equa- 
tions. They can be expressed with respect to the scaled time s = v c t/2, that will 
still be denoted by primes, 

59" (s) - tanh(s)^ / (s) - uj 2 9{s) = (20) 

where 

2 _ 2U"(9 C ) 

{uj can be imaginary). 
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Remark 5. McGehee's equations \18\) are hamiltonian with respect to the sym- 
plectic form a = 2dv A dr 1 / 2 + d^^u) A dO obtained by pullback of the canonical 
form dp x f\dx+dp y /\dy under McGehee transformation. Therefore Morales-Ramis 
applies to this case. 

Theorem 5. Let the Hamiltonian of a system be 

H=^(p 2 x +P 2 y )-U(x,y) 

with U(x,y) > homogeneous of degree —1, defined for all (x,y) / (0,0). Let 
U'{9 C ) = and Let co 2 = 2 g'ff , then if 

uj 2 7^ n{n +1) n being an integer, (21) 

then on a fixed negative energy level the system has no meromorphic integral in a 
neighborhood of the homothetic solution defined by = 9 C . 

Proof. Consider the variational equation (|20p 

d 2 z dz o 

—tt — tanh (t — uj z = 0, 

dt 2 K ' dt 

which is transformed into the differential equation 

d 2 y(t) , . . . , cosh 2 (t) + 4w 2 cosh 2 (t) - 3 , . , . > — - . . 

= 0(t)j/(t), m = W 4cogh 2 (t) ~ , z{t) = y(tW cosh(t). (22) 

Now, by theorem [H the equation (|22ll is ready to be algebrized. The Hamiltonian 
change of variable is r = r(i) = cosh(i), where f = sinh(i), (f) 2 = sinh 2 (t) = 
— 1 + cosh 2 (t) so that 

2 r t — 1 + 4/Ci 

a = -l + T and /= ^ _ ^ . 

The algebraized equation is 

d 2 y(r) r dj/(r) (1 + Auj 2 )t 2 - 3 



dr 2 1-r 2 dr 4t 2 (t 2 -1) 



V(r) = 0, (23) 



and the points 0, 1, —1 and oo are regular singularities. To apply Kovacic's 
algorithm, see appendix A, we use the RLDE 

d 2 V / \ , s 4w 2 r 4 - (6 + 4w 2 )r 2 + 3 , . n 

' - r(r) V , r(r) = A _ 9J _ y \ — 9 , y(r) = ^J== (24) 



dr 2 w " v ' At 2 {t- l) 2 (r + l) 2 ' yK ' 
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T 



with because with u = the differential equation can be solved easily. We 

can see that T = {0, —1, 1, oo} and that the equation (|24p could fall in any of four 
cases of Kovacic's algorithm, now expanding r(r) in partial fractions we have that 

3 8uj 2 - 3 3 3 -8a; 2 3 

r[T) ~ ~16(r- l) 2 + 16(r- 1) ~ 16(r + l) 2 + 16(r + 1) + IT 2 ' 

We start analyzing the case one. The equation ([23]) satisfy the conditions {02,002}, 
because or\ = or_i = 1 = otq = or^ = 2. Due to condition 002, we need the 
Laurent series of r(r) around 00, which corresponds to 

r(r) = u?r 2 + (-\ + At* + (u 2 - f) r 6 + 0(r 8 ) 



obtaining the expressions 

[vU = [Vrl-i = W*\i = [v^]oo = 0, 

at = a+ 1 = l a r=a=1 = I, o± = MvSE. 
By step two, D = Z + and a; 2 has the following possibilities: 

u, 2 = (n + 2)(n + 3), a; 2 = ( 2 "+ 3 X 2 "+ 5 ) , w 2 = (n + l)(n + 2), 
u, 2 = n(n+l), w 2 = (2n+iK2n z i) ) u * = n (n-l), 

discarding to 2 ^ Z because the differential equation has not Liouvillian solutions 
(the monic polynomial P n there is not exists), we take the rest of values for lo 2 
which are equivalents to to 2 = n(n + 1). For each n we can construct u> and by 
step three there exists a monic polynomial of degree n in which each solution of 
the differential equation (|24p is given for all n G Z + . 

Following the case two, we expect to find different values of u> 2 that the pre- 
sented in case one, so that the equation ([241 satisfy the conditions {02,002}, be- 
cause oro = on = or_i = oroo = 2, obtaining the expressions 

E = {-2,2,6}, E x =E-i = {1,2,3}, E^ = {2, 2 - \A + 4w 2 , 2 + y/l + 4w 2 } . 

By step two, Z) = Z + and we obtain again uj 2 = n(n + 1), so that we discard the 
case two. 

Finally, following the case 3, we expect to find different values of uj 2 that the 
presented in case one, but again appear the expression vT+ Auj 2 , which replaced in 
E c and E^ give us again uj 2 = n(n + 1). This means that the differential equation 
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Sorry, the graphic is not available. 

Figure 2: Rectangular 4-body problem 

(|24p is contained in the Borel group when u> 2 = n(n + 1), and it is SL(2, C) when 
lo 2 ^ n(n + 1). Therefore, by remark [21 the Galois group is virtually abelian for 
lo 2 = n(n + 1) and unsolvable for uj 2 ^ n(n + 1). □ 
An alternative proof based on Kimura's approach in given in Appendix A. 2. 

Remark 6. Yoshida \24f gives sufficient conditions for the non-integrability for 
Hamiltonian systems of two degrees of freedom H = (p\ + p 2 )/2 + V(q%,q2) with 
homogeneous potential V of arbitrary integer degree k. He defines the "integrability 
coefficient" A = Trace(Hess V(pi, C2)) — (k — 1) where (01,02) is a solution of the 
algebraic equation 

(c 1 ,c 2 ) = VV(c 1 ,c 2 ) (25) 

and Trace(Hess V(c±, C2)) is the trace of the Hessian matrix. It is not difficult to 
show that the "integrability coefficient" is related to our parameter uj by lu 2 /2 = 
1 — A. Theorem^ can be considered as equivalent to Yoshida's theorem in the 
particular case k = — 1, since condition Ii25\) can be viewed as the vanishing of the 
gradient of the restriction ofV\gi to the unit circle S 1 , i.e. V'(0) = 0. 

Remark 7. Vigo-Aguiar (cited in JEHj) an d co-workers, have developed systemat- 
ically the formulation of Yoshida's result in polar coordinates and used it to study 
two degrees of freedom polynomial potentials. 

Remark 8. The main difference of theorem |UJ and previous work cited in the 
above remarks, is that we are considering explicitly the variational equations along 
a singular ejection-collision orbit. 

In the following subsections we apply the theory developed so far to some 
examples of few body problems. The main interest is to test the non-integrability 
given by theorem [5] in concrete examples having singularities. For simplicity, the 
equations of motion in the examples are recast in McGehee's form (|18p . 

5.2.1 The rectangular 4 body problem 

Four unit masses are at the vertices of a rectangle with initial conditions (position 
and velocity) symmetrical with respect to the axes in such a way that the rectan- 
gular configuration of the particles is preserved. See Figure EJ Let x, y be the base 
and height of the rectangle with the center of mass at the origin, p x = x, p y = y 
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conjugate momenta. The Hamiltonian is 



rr !/ 2 2x 1 1 1 



^ x y y^TF 

Taking polar-like coordinates x = rcos#, y = rsin# the equations of motion are 
of the type (13]) with 

cos 9 sin ft 

The unique homothetic orbit corresponds to 9 C = 7r/4. A simple computation 
shows that 

2 12^2 

u = = 

1 + 2^/2 

then from theorem ((5j) it follows trivially, 

Theorem 6. The rectangular four body problem is not integrable with meromor- 
phic first integrals. 

5.2.2 The anisotropic Kepler problem 

The hamiltonian of the anisotropic Kepler 

1,2 1 



H = -(P Z X +Py) 



depends on the parameter of anisotropy which can be restricted to \x G [0, 1]. For 
[i = and fi = 1 it is integrable. Using polar coordinates x = rcos8, y = rs'mO 
McGehee's equation are obtained (fT8|) with 

U(B) 



yj cos 2 9 + fi sin 2 9 

Homothetic orbits correspond to minima at 9 = 0, tt and maxima at 9 = ir/2, 3ir/2; 
then for minima 

uj 2 = 2(1 -n.) 

According to the non-integrability theorem ([S]), the anisotropic Kepler problem 
is not integrable with meromorphic integrals except when /i = 1 — £(£ + l)/2. 
This leaves only the integrable cases ft = 0,1. Figure [3] shows the Poincare maps 
associated to the section 9 = for some values of the mass parameter. 
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Sorry, the graphics are not availables. 

Figure 3: Poincare maps for the anisotropic Kepler problem (a) \i — 1; (b) 
fi = 0.9; (c) fi = 0.85; (d) [i = 0. Integrable cases \jl — 0, 1 were calculated 
analitically. The main role of two hyperbolic orbits is evident. 



5.2.3 Two uncoupled Kepler problems 

Consider two uncoupled Kepler problems on the line with Hamiltonian 

2 y x y 

which models, for example, two binaries on the line far apart so the interaction 
between them can be neglected. In one binary the particles have the same mass 
taken as unit, and on the other binary \i represents its total mass. This problem 
is evidently integrable for all values of \i. Using polar coordinates x = rcosO, 
y = rsin#, the potential becomes 

cos t> sin t) 

The unique critical point corresponds to the homothetic orbit 9 C = arctan(^ 1 / 3 ). 
One easily computes 

u? = 6 

which does not depend on [i. Thus the non-integrability test ([5|) fails, accordingly 
since the problem is completely integrable. 



6 Open questions and final remarks 

Kovacic algorithm and Kimura's table give the same non-integrability results in 
the specific examples studied in this paper. We recover Yoshida's non-integrability 
in the case of simple mechanical system with two degrees of freedom an homoge- 
neous potential of degre —1, by first performing McGehee's the blow up. Yoshida's 
approach and ours are not entirely equivalent though, since here we are consid- 
ering specifically an ejection-collision orbit exhibiting singularities in the original 
coordinates. This singularity is substituted by an invariant manifold and the sin- 
gular orbit now connects two singular points in the collision manifold. The general 
setting, as stated by Morales [12], of adding singularities to the original Riemann 
surface is here needed in order to apply the theory. To our knowledge only one 
such example is known where this kind of generality is needed, the Bianchi IX 
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cosmologial model, and has been discussed in great detail by Morales and Ramis 
in [15]. Based in this situation we pose the following open problem: 

Consider a Hamiltonian system on a fixed energy level with with an invari- 
ant submanifold A on its boundary. The flow preserves the natural volume form 
on M/j but not necessarily on M^. Let 7 be a heteroclinic (homoclinic) orbit con- 
necting critical points on A but not completely contained in A. We ask: which are 
the class of integrals that are dismissed by Morales-Ramis theory? For example, 
if such an alleged first integral is continuous up to A then by invariance it has to 
be constant on A (following Abraham-Marsden, we call such integrals extendable), 
this is clearly a strong restriction. If such an integral has poles on A are the critical 
points on A necessarily one of them? Another open question is to investigate how 
is transversality of stable and unstable manifolds along 7 related to the solvability 
of the differential Galois group. As a reference, Yagasaki [23] gives an answer to 
this question in the case of A = {p} a critical point an the extended flow Mh is 
still volume preserving (we just add a critical point). 



A Kimura's Theorem 

The hypergeometric (or Riemann) equation is the more general second order linear 
differential equation over the Riemann sphere with three regular singular singu- 
larities. If we place the singularities at x = 0, 1, 00 it is given by 

d 2 j ( \-a-ol + l- 7 -y \ 
dx 2 \ x x — 1 J dx 

x 2 (x — l) 2 x(x — 1) / 

where (a, a'), (7,7')) Q3, P') are the exponents at the singular points and must 
satisfy the Fuchs relation a + a' + 7 + 7' + /? + = 1. 

Now, we will briefly describe here the theorem of Kimura that gives necessary 
and sufficient conditions for the hypergeometric equation to have integrability. Let 
be A = a — a', fi = (3 — j3' and u = u — v' . 

Theorem 7 (pj)]). The identity component of the Galois group of the hypergeo- 
metric equation \26i) is solvable if and only if, either 



(i) At least one of the four numbers X + p, + v, —A + p, + v, A — fb + u, A + fb — u 
is an odd integer, or 

(ii) The numbers A or —A, jl or —p and or —i> belong (in an arbitrary order) 
to some of the following fifteen families 
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1 


1/2 + / 


1/2 + m 


arbitrary complex number 




2 


1/2 + / 


1/3 + m 


1/3 + q 




3 


2/3 + / 


1/3 + m 


1/3 + g 


I + m + q even 


4 


1/2 + / 


1/3 +771 


1/4 + g 




5 


2/3 + / 


1/4 + 771 


1/4 + g 


I + m + q even 


6 


1/2 + / 


1/3 + m 


1/5 + g 




7 


2/5 + / 


1/3 + 771 


1/3 + g 


I + m + q even 


8 


2/3 + / 


1/5 + m 


1/5 + g 


I + m + q even 


9 


1/2 + / 


2/5 + m 


1/5 + g 


I + m + q even 


10 


3/5 + / 


1/3 + m 


1/5 + g 


I + m + q even 


11 


2/5 + / 


2/5 + m 


2/5 + g 


I + m + q even 


12 


2/3 + / 


1/3 + m 


1/5 + g 


I + m + q even 


13 


4/5 + / 


1/5 + m 


1/5 + g 


I + m + q even 


14 


1/2 + / 


2/5 + m 


1/3 + 5 


I + m + q even 


15 


3/5 + / 


2/5 + m 


1/3 + q 


I + m + q even 



Here n,m,q are integers. 



B Alternative proof of theorem [4 

The change of independent variable 

1 

z = - (cos E + 1) . 
reduces the variational equations (|13p to the rational form 



d?£j /3/2 1/2 \ d£ 



I J dz + 2 V-K-z - 1) (-2- l) 2 
By making the (non unique) choice of constants 



1 



1 



(27) 



£j = 0. (28) 



a 



13 = 0, a 
1 



i, 

2' ' 



1. 



7=-(l + V / T+8^ 



7 



^(l-yr+s^) 



4 v- ■ v - ■ - v j ; > 

the equations reduces to the hyper geometric equation of the form given in 

<*j , ( 



d 2 jj ( I- a- a' 1 



7~7 
1 



dz 



■+ 



+ 



77 



1) 



(29) 
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In order to verify the theorem we define the difference of exponents 

/ 1 



A = a — a = —1/2, jX = (3 — (3 = —1, v = 7 — 7' = -yl + 8«J . 

In order to verify condition (i) of Kimura's theorem, we compute the combinations 

1 



A + /t + i> = -(-3 + ^/1 + 8k 



-A + /t + i> = -(-l + yO + 8^ 



v > 



X-jl + U = -(l + y/l + 8Kj) 

1 



X + fi-u = --(3 + 8^1 + 8^-) 

For any of the above quantities to be an odd integer, then Kj must be of the form 

Kj = (n + l)(2n + 3), (n+l)(2n + l), n(2n+l), n G Z. (30) 

In order to verify condition (ii) observe that the only possibility is that /i fits 
in the column of "arbitrary complex number" and A of the form 1/2 + m, with 
m = — 1 an integer, therefore the parameters Kj, j = 1,2 must satisfy the condition 
y/1 + 8kj = 1/2 + £, or 

«i = ^(€ + l). (31) 

But conditions (|30p are contained in condition (|3ip. to see this take I = 2(n + 
1), 2n + 1, 2n, respectively to recover (|30|) . □ 

Remark 9. H^e recovered condition (14\ ) 



C Alternative proof of theorem [5 

The change of dependent variable 

56{s) = cosh 1/2 (s)y(s) 

reduces the equation to 

y"(s) - -(1 + 4w 2 - 3sech 2 (s))y(s) = 0. 

A further change of independent variable z = sech 2 (s) yields 

dz = — 2 sech 2 (s) tanh(s) ds = —2zyl — z ds. 
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Thus 



d 2 y r d ( r dy 

— = Az\/l - z— z\J\ - z— 
as az \ az 



developing the second derivative 

4z 2 (l - z)y"{z) + ( 4*(1 -z)+ Az 2 Vl^z~ 



-1 



4z 2 (l - z)y"{z) + 2z (2(1 - z) - z) y\z) 



(1 +4uj 2 - 3z)y(z 
(1 +Aoj 2 - 3z)y(z 



4z 2 (l - z)y"{z) + 2z(2 - 3z)y'{z) - -(1 + 4uj 2 - 3z)y(z 



z{l-z)y"{z) + -(2-3z)y'{z) 



16z 



(1 + Aoj 2 - 3z)y(z 



y"(z) + 



1 - p 
zO—z) 



y\z) + 



1« 4 16 



z 2 (l-z) 



Expanding in partial fractions we finally get 



//'(--) + I - + 



1/2 



z z — 1 



, 1/16-kj 2 /4 -l/8 + w 2 /4\ 



z(z-l) 



which is a Riemman equation. Comparing with (|29p a convenient choice of pa- 
rameters is 7' = and 



a 



--y/l + 4w 2 , a' = -Vl + 4w 2 



The exponent differences are 



A 

A 



a — a 

1 
1 

2 



-^Vl + 4w 2 



Condition (i) of Kimura's theorem is satisfied whenever any of the four combina- 
tions indicated there is an odd integer; thus 



uJ- 



n(2n — 1), n(2n + 1), n being an integer 



(33) 



To verify condition (ii) of Kimura's table, notice that p, = — 1 is not of any of the 
forms of the columns except for the first case: We can take fi as an "arbitrary 
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complex number" and v = 1/2 of the form 1/2 + m, with m = 0; thus in order to 
fit the first case A must be of the form 1/2 + 1, I an integer, that is 




this yields the condition 



£{£ + !) 



(34) 



Now observe that condition ([33]) is contained in condition (|34f) by taking I = 2n — l 
or I = 2n. 



Remark 10. We recovered condition [2l\) . 
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